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1 Introduction 



Recently there has been considerable interest in the novel suggestion that we live in a 
Universe that possesses more than four dimensions. The standard model fields are assumed 
to be confined to a hyper-surface (or 3-brane) embedded in this higher dimensional space, in 
contrast to the gravitational fields which propagate through the whole of spacetime [HO]. 
In order for this to be a phenomenologically relevant model of our universe, standard four- 
dimensional gravity must be recovered on our brane. There are various ways to do this, 
the most obvious being to assume that the extra dimensions transverse to our brane are 
compact. In this case gravity can be recovered on scales larger than the size of the extra 
dimensions [5-7]. This is different from earlier proposals since the restrictions on the size 
of the extra dimensions from particle physics experiments no longer apply, as the standard 
model fields are confined to the brane. The extra dimensions only have to be smaller 
than the scale on which gravity experiments have probed, of order 0.1mm at the time of 
writing. Another way to recover four-dimensional gravity at large distances is to embed a 
positive tension 3-brane into an AdSs bulk [9,10]. In this scenario four- dimensional gravity 
is obtained at scales larger than the AdS radius. Randall and Sundrum showed that this 
could produce sensible gravity even if the extra dimension was not compact. 

Several interesting aspects of the above extra dimensional scenarios have since been 
investigated, and compared with the standard four-dimensional case. The cosmology of a 3- 
brane in a five- dimensional bulk was studied and its Friedmann equation derived and shown 
to contain several extra terms [11-15]. Perturbations to this homogeneous case have been 
examined [16-19] as have some infiationary models [20,21], as well as alternative solutions 
to the fiatness and horizon problems [22]. The behaviour of phase transitions, topological 
defects and baryogenesis in a brane world setting were also considered [23] . More recently 
it was shown how to embed the Randall-Sundrum models within supergravities [24-26] 
and then within string theory compactifications as in [27-30]. The seemingly arbitrary 
feature of having an AdS bulk spacetime is actually well motivated as it is found as a 
super symmetric vacuum to supergravity theories, inspiring several more recent brane world 
models [31-37]. For recent, comprehensive reviews of the subject see [38,39]. 

The extra term C/r"^, in the Friedmann equation known as the dark radiation term, 
could provide a possible test for the above extra-dimensional scenarios. Although many 
authors have assumed that C is constant, this will not be the case if the brane emits 
gravitons at early times. The evolution of C has been studied [40-42] for symmetric branes. 
However, in this paper we examine the dynamics of a non-Z2 symmetric graviton emitting 
brane [43] and determine what consequences the lack of symmetry has on the evolution of 
C. We begin in section 2 by summarising the main aspects of non-Z2 branes, including the 
new Friedmann equation and the global nature of the five- dimensional spacetimes either 
side of the brane. Then in section 3 we derive the equations that govern the evolution of 
the energy density p, the dark radiation parameter C and the non-Z2 symmetric parameter 
F, using a Vaidya spacetime approximation. The global nature of these five-dimensional 
(y4(iS')-Vaidya spacetimes is then discussed along with the appropriate conformal diagrams. 
We then analytically solve the system of equations in the high/low energy limits in section 4, 
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before presenting a full numerical treatment in section 5. We examine the full range of 
initial conditions for C and F, and include several interesting cases that have so far been 
ignored. Constraints on the size of the initial dark radiation and non-Z2 symmetric terms 
in this scenario are derived and are notably independent of M5, the five dimensional Plank 
mass. Our conclusions are given in section 6. 

2 Non-Z2 Symmetric Branes Moving Through Five- 
dimensional Spacetimes 

In this section we summarise the relevant aspects of a non-Z2 brane in an AdS bulk. We 
discuss the derivation of the non-Z2 Friedmann equation and also the global nature of the 
five-dimensional Schwarzschild-Anti-deSitter spacetime. 

2.1 The Non-Z2 Symmetric Friedmann Equation 

Starting with the five- dimensional Einstein equations, 

Gab + ^9ab = i^^Tab, (1) 

where k is related to the fundamental five- dimensional Plank mass M5 by k = 1/M|, one 
can use Birkhoff 's theorem to show that the bulk metric must take the form, 

ds^ = -f^{r)dT^ + r^dJ^l + J^^d^^- (2) 

Here dJ^l represents the maximally symmetric three-dimensional metric with A; = — 1,0, 1 
parameterising the spatial curvature. The ± in equation 2 represents the fact that the 
bulk spacetime can be different either side of the brane. The function /^(r) is found by 
substituting the metric (2) into the Einstein equations (1) and is given by [15], 

/±(r) = /iV^ + k - (3) 

fi is the inverse curvature radius, related to the five-dimensional cosmological constant by 
/i = ■y/— A/6, and is the bulk black hole mass parameter to the right /left of the brane 
respectively. We see that the bulk solution on either side of the brane can be either an 
AdS black hole (C^ > 0), an AdS naked singularity (C^ < 0) or pure AdS (C^ = 0). 

If a brane with stress-energy-momentum tensor T^b = 5(r)diag(— p — a, P — a, P — 
a,P — a, 0) exists between two such spacetimes, the Israel junction conditions [44-47] can 
be used to show that the Friedmann equation governing the expansion of the brane is given 
by [48], 

~ ~ 18^ ^ 36^ ^2 + ^4 + (p + ^)2^8- ( ) 
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Here a is the bare brane tension, p is the physical energy density of the brane, and we have 
imposed the Randall- Sundrum tuning conditions such that + A/6 = 0, and hence 

that k'^(T^/36 = /i^. The Weyl (or dark radiation) parameter C and the non-Z2 parameter 
F are dimensionless, and are related to the average and difference of the black hole masses 
either side of the brane via, 

C+ + C+ - C~ 

^ - 

The implications of the extra terms in the Friedmann equation (4) have been thoroughly 
investigated [11,12,48,49]. However, as will be discussed in the following section, if gravi- 
tons are emitted by the brane into the bulk at early times then the parameters C and F will 
no longer be constant and the early cosmological evolution of the brane could be greatly 
affected. Note that there is an alternative method of breaking the Z2 symmetry across 
the brane; one can have different cosmological constants either side of the brane leading 
to a term oc 1/r^ in the Friedmann equation. Although interesting, this type of non-Z2 
symmetry is not directly affected by graviton emission and so will not be considered here. 



2.2 Global Structure of Sch-AdS 

Here we summarise some aspects of the global structure of the aforementioned Sch-AdS 
spacetime and the brane 's trajectory through it. For clarity we refer only to the spacetime 
on the right hand side of the brane, which has black hole mass . Such concepts will be 
of use when in the next section we generalise the analysis to Vaidya spacetimes. In what 
follows, we assume the case where /c = and is positive. Other cases can be argued 
similarly. 

One should note that the coordinates T and r in equations (2) and (3) do not cover the 
whole of the Sch-AdS spacetime, and are valid only for r > = {C^Y^"^/ p^^^. To continue 
beyond this horizon we first note that the metric given by equation (2) is conformally 
equivalent (after rescaling) to, 

ds' = - (^) rfT^ + dr- + r^d^l (6) 

where the horizon is now at r = 1. We then map the radial coordinate from 1 < r < 00 to 
— (X)<i?<7r/4 according to the transformation. 



, r'^dr 1 . , ^ 1 , 

R = I = -tan-i(r) + -log 



r — 1 

r + 1 



(7) 



Then after changing to the null coordinates given hju = t — R and v = t + R the metric 
looks like, 

ds' = - ( ^—^j dudv + r^dT.1. (8) 
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Finally changing to the coordinates U and V given by, 



U 



V 



— exp(— 2m), 
exp(2?;). 



(9) 
(10) 



results in the metric. 




dUdV 
-4UV 



) 



+ r'dj:l 



(11) 



Note that the point (1, 0) in the (r, t) plane has now been mapped onto (0, 0) in the {U, V) 
plane and hence that the future and past horizons at r = 1 of the black hole have been 
mapped onto the positive V and negative U axes respectively. To see that the zero at r = 1 
has been removed note that. 



which behaves as (r — 1) for r ^ 1. Using the above coordinate transformation one can 



Figure 1: Conformal diagram of the extended Schwarzschild Anti-DeSitter spacetime. The 
thick curved line represents the trajectory of a cosmologically realistic brane, while EH 
labels the event horizon. Note that here, and in all subsequent figures only the spacetime 
to the right of the brane exists. 

deduce that the conformal diagram of five- dimensional Sch- AdS is given by figure 1. The 
left and right vertical boundaries represent the spatial infinities, while the top and bottom 
boundaries correspond to the black and white hole singularities respectively. Although 



UV 




(12) 



r = 




r = 
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the Friedmann equation (4) is derived from the metric (2), it is still valid in all regions 
of the Sch-AdS spacetime, and hence we can use it to examine the trajectory of a brane 
moving through such a spacetime. A cosmologically realistic brane that starts from an 
initial singularity r = with C"*" > 0, and then expands forever will have the trajectory 
shown in figure 1. It begins at the white hole singularity, passes through the event horizon 
at r = 1 and continues off to infinity. 

We can ask how the situation would change if the brane is emitting gravitons toward 
the bulk, and more specifically how this would effect the black/white hole mass parameter 
C~^. In the next section we use a five-dimensional AdS generalisation of the ingoing Vaidya 
spacetime to model the effect of graviton emission by the brane, and hence derive the 
evolution equations for and . We then go on to discuss the extension of Vaidya 
spacetimes, and present suitable conformal diagrams to explain conceptually the results 
given in sections 4 and 5. 



3 Graviton Emission and Vaidya Spacetimes 

In the previous section we discussed a brane moving through two different spacetimes 
possessing only a cosmo logical constant. We now go on to generalise the situation presented 
by Langlois et.al. [41], where the brane emits gravitons at early times, to the above non-Z2 
scenario [43]. 

In order to model this situation we use the five-dimensional generalisation of the in- 
coming^ Vaidya metric, 

ds"^ = -f^{r,v)dv'^ + 2drdv + r'^d^l, (13) 
where we once again assume a different spacetime either side of the brane and that k = 0, 



fHr,v) = A'-^. (14) 



Here i; is a null ingoing coordinate. Note that now the Weyl parameters C^{v) no longer 
represent purely the mass of the black hole in each spacetime, but instead represent the 
energy contained inside the radius r for a given v, and hence describes an inward flow of 
radiation. If did not depend on v, then the metric (13) would just be a rewriting of 
Schwarzschild-AdS, given by equation (2), as can be seen by the coordinate transformation 
V = T + J dr/f{r). The brane moves through both such spacetimes with trajectories given 
by v^{t) and r(t), where t is the time experienced by a brane based observer. Normalisation 
of the brane's velocity vector in each spacetime u"^^ = {u^, r, 0, 0, 0} requires that. 



r 



r2 



f 



± 



^In the bulk the brane resides at the largest radius of the spacetime, hence any gravitons emitted will 
travel inward. 
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which we will need later. 

The incoming Vaidya metric is a solution to Einstein's equations with a bulk energy- 
momentum-tensor corresponding to null ingoing radiation, 

T^B = i^kik% (16) 

where k^"^ are null ingoing vectors in each spacetime. By inserting the metric (13) and the 
stress-energy-tensor (16) into Einstein's equations we find that they will be solved provided 
the Weyl parameters satisfy, 

dC^v) _ 2^ 

The appropriate normalisation of k^"^ is given by k^u"^^ = 1, which implies that the only 
non-zero component is k^^ = k^ = Also ip is the flux of gravitons leaving a radiation 

dominated brane as seen by a brane observer, and has been shown by Langlois et. al. [41] 
to be. 



a 



^ = (18) 

where a is a dimensionless constant. Note that here we implicitly assume that gravitons 
are produced in equal amounts either to the left or right of the brane. Particle interactions 
that produce these gravitons will not feel the effect of different bulk masses either side of 
the brane. However we will go on to show that the difference in the brane's trajectory in 
each spacetime alters the effect of the graviton emission so that it appears to be \10n-Z2 
symmetric. 

Combining equations (15), (17) and (18), we get an expression for the evolution of 
either Weyl parameter due to graviton emission, 

dC^ an^ 3^2^ 



dt 



-r"p\^/r^ + -r). (19) 



In order to determine the behaviour of a brane positioned between two such Vaidya space- 
times, we need to impose the Israel junction conditions [44-47] across the brane. These 
conditions relate the jump in the extrinsic curvature of the brane K^b = c^Bi where 

as usual is the unit vector normal to the brane and Hab = Qab — nAfiB the induced 
metric, to the brane's energy momentum tensor tab = diag(— p — a, P — a, P — a, P — a,0). 
This can be written in the covariant form: 

[Kab] = t^\rAB - ^rhAB). (20) 

The ordinary spatial components of these Junction Conditions show that the Friedmann 
equation for the brane is of the same form as that of equation (4). However, now both the 
Weyl parameter C and the \1011-Z2 parameter F are now functions of v and hence t. Using 
the fact that H = r /r, we can use equations (4) and (3) to rewrite equation (19) as. 



di 18 V 6 2^2 (p + a)r^ 
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It can immediately be seen that the process of graviton emission will lead to a reduction of 
the Z2 symmetry since if C"^ is greater/less than C~ then dC^ /dt will be less/greater than 
dC~ /dt, and hence the non-Z2 parameter F will decrease. To analyse this in more detail 
we take the sum and difference of equations (21) to obtain expressions for the evolution of 
the Weyl parameter C and the non-Z2 symmetry parameter F defined in equations (5), 



dC 


A '•'2 A / 

an p r ( 


di 


18/i2 V 


dF 




di ~ 


6(p + a) 



y(p + ^)-^), (22) 

F. (23) 

These equations will be examined both analytically and numerically in sections 4 and 5 
for various cases of interest. 

The time and fifth dimension part of the Junction Conditions, given by equation (20), 
also yield another relation representing the non-conservation of energy on the brane, 

'p + 3-(p + P) = -2^. (24) 
r 

The factor of 2 on the right hand side is due to the fact that the brane is emitting a 
flux ip of gravitons into each spacetime either side of the brane. We have hence derived 
the four 'master' equations that describe brane graviton emission into non-Z2 symmetric 
spacetimes, given by (4), (22), (23) and (24). In the next section we investigate the global 
properties of Vaidya spacetimes in order to understand the results presented in sections 4 
and 5. 



3.1 Extension of the 5D Vaidya- AdS Spacetime 

We have derived the relevant equations that describe both the motion of a graviton emit- 
ting non-Z2 symmetric brane and the evolution of the mass parameters either side of the 
brane. However, one should note that the coordinates [v, r) do not cover the whole of 
the spacetime. Therefore we now show how to extend the Vaidya coordinates (see [50,51] 
for four dimensional examples), in order to gain a conceptual understanding of certain 
cases of interest, including when the brane emerges from the white hole. Previous authors 
have ignored these solutions, mainly by assumption [41,43,52]. We mainly concentrate 
on the situation where the BH mass parameter C'^{v), is positive for all time; cases where 
C~^{v) < or where C~^{v) changes sign can again be argued similarly [51]. 

The ingoing five-dimensional Vaidya metric (13) is conformally equivalent (after rescal- 
ing r — >• r/yU and v ^ v/ p) to: 

ds"" = -(^^- dv^ + 2dvdr + r^d-Ll (25) 

where rf^{v) = pi?C^{v). Examining the Friedmann equation (4), it can be seen that if 
the BH mass parameter is initially greater than zero {C^{v — > —00) > 0), the brane 
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will originate from the singularity inside the white hole (r = 0) and emerge from it at a 
later time, similar to the non-emitting case. The above Vaidya coordinates only cover the 
interior region of the black hole (along with one of the exterior regions) and hence need 
to be extended. The conformal diagram of the region covered by the {v, r) coordinates is 

r = 




Figure 2: Conformal diagram of the unextended ingoing AdS Vaidya spacetime 

shown in figure 2, along with a possible trajectory of the brane (r = a(t)), the direction of 
graviton emission (the dashed arrows) and the various horizons. EH represents the normal 
event horizon of the black hole, whether PEH represents the past event horizon appearing 
at r = rh{v — ^ — oo). It can be seen that outgoing radial null geodesies emerge from 
the PEH at finite values of their afiine parameter and are hence incomplete: extension is 
therefore required. 

Figure 2 also shows one of the main features of non-stationary spherically symmetric 
spacetimes, which is the decoupling of the apparent horizon (AH) from the event hori- 
zon [50]. The apparent horizon is the hypersurface separating the regions with or without 
closed trapped surfaces (see [53] for more details), here given by r = rh{v), and will be 
matched to a corresponding AH in the extended spacetime. 

In general, extending Vaidya spacetimes is non-trivial as is discussed by [50]. One 
usually does not know the form of the 'mass function' C~^{v) (and therefore r^iv)) in the 
extended part of the spacetime or indeed even whether radiation is ingoing or outgoing 
there. This leads to many possible choices of extension. However, in the case whereby a 
brane emits the radiation into the bulk, we do in principle know the general form of C+(f), 
and hence we can adapt Israel's extension of four-dimensional Vaidya spacetimes [54] to 
the five-dimensional AdS Vaidya metric given in equation (25). 

Changing coordinates from (t',r) to (f/, V) by the transformation: 

dx 

r = AU + rh, V = —- (26) 
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where the function A(y) is chosen to be: 



A{V) = I 2rh{x)dx, (27) 
Jo 

the metric given by equation (25) then becomes, 

= ^\w + rj"^ +1^) + + (AU + r,fdY.l (28) 

Note that the PEH has been mapped from (r = r^^v = — oo) to (—00 < U < 00, V = 0), 
and AH from (r = r/i(f ), —00 < v < 00) to (t/ = 0, < V < 00). The extended part of the 
spacetime corresponds to —00 < V < 0. Some conditions must be imposed upon the mass 
function rh{V) to ensure the finiteness and continuity of both the metric and the energy 
momentum tensor [54] (which we know to still describe radiation leaving the brane), which 
are that rh{V) must be a function such that [50], 



MV = 0) ^ 0, ^"'^ ^ 



dV A 



v=o 



< 00 ^ ^{V = 0) = 0. (29) 

(Jj V 



We have assumed the first condition to be true previously. As for the second condition, 
demanding that the full metric satisfies the energy conditions [54] implies that, 

drh r > for y > 

dV \ < for y < ^ ^ 

and hence the second condition is satisfied purely by demanding that the metric is physi- 
cally realistic. This has some interesting implications: we can see that when the brane is in 
the extended part (V < 0) of the spacetime, any graviton emissions will actually decrease 
the size of the mass function rh{V), and hence decrease C~^(y). The conformal diagram 
for the extended spacetime is given by figure 3, which also shows the the direction of the 
emitted gravitons and the trajectory of the brane. This trajectory is a possible solution 
r{t) to the equations (4), (22), (23) and (24) which are still valid in the extended part of 
the spacetime. One can see that the V < region is actually a region of outgoing Vaidya 
spacetime with the metric. 



ds' = - ( f2 _ !lM ) dv^ - 2dvdf + r^dEl (31) 



One should also note that the coordinates {U, V) used to extend the spacetime are only 
valid provided AU < \/2rh, which means they cover the interior of both the black and 
white holes, but do not extend to spacelike infinity (r = 00). 

Figure 3 shows that any gravitons emitted before the brane crosses V = (equivalent to 
r = Tfi) are sent in an outward null direction and hence leave the white hole, subsequently 
escaping to infinity, and will therefore cause the mass parameter C'^ to decrease. This is in 
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r = 



Figure 3: Conformal diagram of the extended ingoing AdS Vaidya spacetime 

agreement with equation (21) which gives that dC'^ /dt will be greater/less than zero if r is 
greater/less than (and a similar result for C~). This situation where the mass parameter 
initially decreases will occur at early times in the brane's cosmological evolution, when the 
brane will be at high temperature. This raises the possibility that graviton emission at 
early times could help to decrease (and/or C~ and therefore also C), and hence decrease 
the amount of dark radiation below the bounds set by the abundance of light elements 
at nucleosynthesis. Figure 3 shows the case where the mass function is symmetric either 
side of \^ = 0, but for an expanding, cooling brane thermally emitting gravitons the 
situation would be very asymmetric and could lead to the final mass of the black hole 
being substantially smaller than that of the white hole. 

The case where the black hole mass parameter is initially negative, before changing to 
positive is shown in figure 4. Initially, < and there is a naked singularity at r = 0. 
However, provided the brane emits enough gravitons, will change sign and a black hole 
will be formed [51]. In this situation, equation (21) shows that always increases. One 
must of course note that a naked singularity in the bulk is phenomenologically problematic 
for a realistic brane world model. We can now use figures 3 and 4 to understand the general 
case. As has been discussed above, regardless of the initial values of and C~, the non- 
Z2 parameter F will always decrease. The behaviour of the dark radiation parameter C 
is as follows: if < C~ < C"*" (the spacetimes either side of the brane both correspond to 
figure 3), then C will initially decrease until some time where C~ < /x^r^ < C~^, from which 
point C will increase forever, asymptotically approaching a constant value as t — * 00 and 
the brane cools. If C~ < < 0, (the spacetimes either side of the brane both correspond 
to figure 4), then C will increase for all time, asymptotically approaching a constant as 
before. If C~ < < C~^, (the spacetime on the right/left side of the brane corresponds to 
figure 3/4), then initially C~ will increase, C"*" will decrease and the exact behaviour of C 
depends on the relative sizes of and C~ . In each of the above cases the final magnitude 
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Figure 4: Conformal diagram of a naked singularity changing to a black hole 



(and sign) of C is as yet undetermined so in the next two sections we give analytic and 
numerical solutions that can be understood in terms of the conceptual arguments presented 
in this section. 

Throughout this paper we have assumed that the gravitons are emitted in a perpen- 
dicular direction from the brane and hence that a Vaidya spacetime is a realistic model 
of the situation. Examining the case where gravitons are emitted in all directions is far 
more complicated [55], and is left to future work. However, we can use the above con- 
formal diagrams to gain some understanding of the effect of such emissions: in figure 3, 
non-perpendicular gravitons would follow timelike geodesies leaving the brane. Therefore, 
when the brane is inside the white hole, not all of the gravitons that it emits during this 
time will escape to infinity, and hence the mass parameter would not decrease as much as 
is described both above and in the following two sections. 

4 Equations for the Evolution of a Graviton Emitting 
U0U-Z2 Brane 

We now go on to discuss analytic solutions valid in the high/low energy regimes and when 
the C and F terms in the Friedmann equation are assumed to be subdominant. More 
general situations will be examined numerically in the next section. 
First it is useful to define the following dimensionless parameters, 

P - H 

p = t = fit, H = —. (32) 

a n 
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The equations (4), (22), (23) and (24) can then be written in a somewhat simpler dimen- 
sionless form, 



= p" + ^p + ^^ + i^TW'' ^^^^ 

I = -AHp-ap^ (34) 

^ = 2ar'p\p + l-H), (35) 

dF ap^ 



dt p + 1 



F. (36) 



Unfortunately, finding a general analytical solution to these equations, such as that found 
by Leeper et. al. [42] for the Z2 symmetric case, is very difficult and hence we restrict our 
analytical investigation to examine the early and late time limits only. 



4.1 The Early Time/High Energy Limit 

Assuming a high energy limit such that p ^ 1, and also that the C and F terms in the 
Friedmann equation are sub-dominant, we can solve the evolution equations approximately. 
Under these assumptions the Friedmann equation (33) takes the much simpler form of 
~ p^. Combining this with equation (34) and integrating we obtain. 

Equation (34) also shows that p = ■y/r^~^°', where 7 is a constant, and hence that, 

r(t) = [7(4 + a)t]3T^, (38) 

as opposed to the usual r ~ t^^^ in early brane world cosmology without graviton emission. 
To find the time dependence of the mm-Z^ parameter F we use the fact that p ^ 1 and 
hence equation (36) becomes dF/dt ~ —apF, which when integrated using equation (34) 
gives, 

F = Ap^ = (39) 

[(4 + a)t]3+^ ^ ' 

where A is a constant set by the initial conditions. It can be seen that although, as 
expected, F decreases in this limit; it does so at a much slower rate compared to p. What 
is important here is the behaviour of the F term in the Friedmann equation (33), and we 
can use the above results to show that, 

F^ 

(40) 



(p + l)2r8 78/(4+") • 
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Since the F term behaves as a constant this raises the possibihty of it dominating the 
Friedmann equation for a period, similar to the case without graviton emission [48]. To 
find the time dependence of C{t) we need to use the following approximation of H, 

Inserting this into equation (35) and using equation (40) we obtain, 

- -V(l-^-#V (42, 



dt 7^/(^+") H 

The C{t)/r^ term on the right hand side should not, in general, be neglected as was done 
in [41]. Using the above expressions for p{t) and r(t), this becomes, 

dCitl ^ ^ _ aCjt) 
dt (4 + a)t' ^ ' 

where i? is a constant given by. 



4 , 4 

/4,+a 



B = a\- (44) 

(4 + a)4+<. 



Equation (43) has the solution. 
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C{t) = B[t^ + (45) 

where to is an integration constant. Again, the important fact here is the behaviour of the 
C{t)/r^ term in the Friedmann equation, compared to the dominant term. This can be 
written as, 

^ = -^(l-AV7^/^+")(l+to(4 + a)p). (46) 

So in general this term is not initially constant as assumed by [41], but instead has a p 
dependence, which could lead to it dominating the Friedmann equation at a later time, as 
will be discussed in the next section^. This dependence agrees with the exact Z2 symmetric 
solutions found in [42]. Note the effect here of the non-Z2 symmetric term: aside from 
acting like a cosmological constant in the Friedmann equation it has a significant effect on 
the evolution of the Weyl parameter C. For example, if initially C = and ii A > 7^/^+" 
then C will behave as C ~ —r^ i.e. will grow more negative. This corresponds to the case 
discussed in section 3.1, where C~ < < and there is initially a naked singularity /black 
hole on the left / right side of the brane. If C remains negative this could result in if = 
which would lead to the universe collapsing, and has to be investigated further. 



^The derived form of both the C and F terms in the Friedmann equation vahdate the assumption that 
they are sub-dominant in this regime 
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4.2 The Late Time/Low Energy Limit 



For the low energy/late time regime, we assume both that p ^ 1 and again that the C 
and F terms in the Friedmann equation are sub-dominant such that now ~ 2p. Using 
this, equation (34) can now be integrated to give. 



1 



Pi 



■1/2 



+ 2V2{t - U 



2 ' 



(47) 



Note 



where p(tj) = Pi, tj being the time the low energy approximation becomes valid, 
the usual p ~ dependence of standard four- dimensional cosmology. Equation (34) also 
shows that r = (7/p)^''^, and inserting this into equation (36) and integrating, one obtains 
the solution for the non-Z2 parameter in this late time regime in terms of p. 



F{P) 



Fjexp 



-f?." (, _ (p: 



6^2 



Pi 



(4J 



where F{pi) = Fi. Since the energy density at which the low energy approximations become 
valid is given by pi ~ 2, we see that in this regime F roughly decreases by a factor of e~"/^ 
which is small as expected since the universe is now cooler and hence emits less gravitons. 
Replacing the above results for p(t), r{t) and F{t) into equation (35) and integrating; we 
get the solution for the behaviour of C in the late time limit which we write in terms of p 
for convenience, 



C{P) 



1-1^ 

pt 



1/2 



1 



Pi. 



(49) 



where C{ti) = Ci. Depending upon the assumed value of pj, we can see that in this limit C 
roughly increases by 07, a small amount as expected. Note that in this regime the \10n-Z2 
symmetry has no effect on the evolution of C, mainly due to being rapidly damped away as 
r increases. To analyse the effects of larger F and C terms, and to examine what happens 
between the regimes, we now turn to numerical methods. 



5 Numerical Analysis 

Previously we have studied some aspects of non-Z2 symmetric graviton-emitting branes in 
certain limits, but in order to gain a full understanding of the effects of non-Z2 symmetry 
in this scenario we solved the equations (33), (34), (35) and (36) numerically for a range 
of values of both C and F and the results are shown in figures 5 and 6. 

Figure 5(a) shows the time dependence of the four contributions to the Friedmann 
equation: p^, 2p, C/r^ and F^/{{p + l)^r^), on a log scale, with the initial values (labelled 
by a subscript z) taken to be 2pi = 10000, Ci/rf = 500 and Fi/{{p, + l)rf) = 0.05. The 
initial time ti was chosen to be in agreement with equation (37) such that tj = ((4 + a;)pj)~^ 
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in order for t = to be the time of the initial singularity. The first thing to notice is the 
behaviour of each of the terms as t 0. One can see that the term dominates as 
expected, the F term tends to a constant in agreement with equation (40), and the C term 
is proportional to p (as predicted by equation (46)) and does not necessarily tend to a 
constant as assumed by [41]. As t ^ oo, the brane cools and graviton emission is hence 
reduced. Therefore, as in the non-emitting case, the F term is proportional to p^, and C/r^ 
is proportional to the dominant 2p term. Note that the ratio of dark radiation to normal 
radiation defined by, 

- f^. (50) 

increases only during the transition between these two limits, hence numerical analysis is 
definitely required to determine the final dark radiation content of the brane world model. 

The evolution of C, 2F, and C~ is shown in figure 5(b), using the same initial 
conditions as in the previous figure. Note the behaviour of C+: initially it decreases, as 
the brane has yet to emerge from the white hole on its -|-ve side. This corresponds to the 
extended part of the Vaidya spacetime as shown in figure 3, where any graviton emissions 
actually leave and hence reduce the size of the white hole. Once the brane passes the 
horizon on its +ve side, then begins to increase. C~ increases from zero in a similar 
manner to the limited cases discussed by [43]. Figure 5(c) shows a similar situation except 
the initial value of C (and hence C~) are both negative. The initial conditions are the same 
as the previous diagrams apart from now Ci/rf = —400. One can see that C~ changes 
from negative to positive. This means that on the — ve side of the brane there is initially a 
naked singularity (C~ < 0) which later turns into a black hole {C~ > 0). This corresponds 
to the conformal diagram given in figure 4, where the incoming gravitons always cause the 
mass parameter to increase. In general, all possible behaviour of F and C can be explained 
conceptually by figures 3 and 4. 

The evolution of C(t) for different initial values of the F term (ranging from to 
1.5) are shown in figure 5(d). It can be seen that at early times C(t) decreases with the 
C{t) ~ behaviour given by equation (45), while at late times when there is less 

graviton emission C{t) tends to a constant; the specific value of which is lower for larger 
initial F. One can see that the conclusions made by [43], that a non-zero F can help reduce 
the final value of the dark radiation term and hence satisfy the nucleosynthesis bounds, 
are correct for small initial values of F. However, what has not been realised previously is 
that if F is too large then, as is shown in figure 5(d), the final value of C{t) can be negative. 
This raises the possibility of not only violating the lower nucleosynthesis bound on C, but 
also of leading to H = 0, at which time the universe begins to collapse. An example of this 
is given in figure 6(a) which has initial values Ci/rf = 5000 and Fi/{{pi + l)rf) = 22, and 
shows the magnitude of the contributions to the Friedmann equation on a log scale. One 
can see that at a time given by log(t) ~ —6.1 the C term changes from positive to negative, 
and its magnitude subsequently becomes larger than the 2p term. Since at late times the 
F term decreases much faster than the C term, the latter term will dominate. This leads 
to if = at a time given by log(t) ~ 0.56, and hence the universe collapses. Note that 
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previous analysis of non-Z2 symmetry in the non-emitting case found that a large F term 
was phenomenologically acceptable (at least in terms of the expansion rate) as this term is 
rapidly damped at late times. However, now we can see that this is not the case if the brane 
emits gravitons, as a large F term leads to the universe collapsing. Figure 6(b) shows the 
case where the non-Z2 symmetry actually helps satisfy the nucleosynthesis bounds. Initially 
Ci/rf = 50 and Fi/{{pi + l)rf) = 2 and here the dark radiation term changes from positive 
to negative and then back to positive again, and the final ratio e/ is given by, 

ef = e{t ^oo) < 0.037, (51) 

which is well inside the current constraints. For convenience we now define the value of 
the F term as t — as, 



(T) 



(52) 



(p + l)r4 

(T) 

Note that both F^ and e/ have well defined constant values. 

In order to understand the full effect of non-Z2 symmetry in this scenario, we graphed 
the dependence of the final ratio c/ on the initial size of the F term given by F^ for 
values of the emission parameter a = 0.02, 0.2, 0.4, and 1.0 as is shown in figure 6(c). In 

(T) 

all cases the initial size of the dark radiation parameter C was zero. When is small, 
the larger a is the larger ef will be, as expected. However, one can see that as discussed 

(T) 

above the larger the magnitude of F^ , the smaller the final ratio ef. Ife/<— 1 this will 
inevitably lead to = at some point in the late time regime (when only the 2p and C 
terms can dominate) and hence cause the universe to collapse. Note how a larger a leads 

(T) 

to a much smaller range of allowed values for F^ that ensure that ultimately e/ > —1. 

Figure 6(d) again shows a graph of the dependence of the final ratio e/ on the initial 
size of the F term: F- . However, now the emission parameter is set to a = 0.02 (as 
calculated by [41]) and the initial ratio = {Ci/rf)/(2pi) takes the values -0.25, 0, 0.25 
and 0.5. The two horizontal lines represent the upper and lower nucleosynthesis bounds 
on the ratio of dark to normal radiation [56] given by, 

-0.41 <ef < 0.105. (53) 

We can therefore obtain a combined restriction on the initial size of the F and C terms. 
The results in figure 6(d) can be accurately approximated by the expression: 

ifPi 



+ + ei. (54) 

205 4 + a ^ ' 

The first term on the right hand side represents the effect of the non-Z2 symmetry, while 
the second term is the approximate increase in the ratio of dark to normal radiation for 
a Z2 symmetric brane as calculated by [41]. Combining equations (53) and (54) gives the 
following restriction on the magnitude of the initial F term in terms of the initial ratio ej 
due to the nucleosynthesis constraints, 

-20.5 + 205 < < 85.1 + 205 e^, (55) 
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where we have kept all numerical terms to three significant figures. This shows that if 
the initial dark radiation term is more than 10% of the normal radiation content of the 
universe, such that > 0.1, then < |-Fj | < 106 and hence the brane has to be non-Z2 
symmetric; otherwise the nucleosynthesis constraints would not be satisfied. If on the other 
hand < —0.415 then the right hand side of equation (55) is less than zero and regardless 
of the size of the F term, the nucleosynthesis constraints will be violated. 

6 Conclusions 

We have analysed the dynamics of a non-Z2 symmetric brane that emits gravitons at 
early times, using an AdS-Vaidja. spacetime approximation. The equations governing 
the evolution of the energy density p, the scale factor a, the dark radiation parameter C 
and the non-Z2 symmetry parameter F were derived and analysed. We then discussed 
the structure of the five- dimensional AdS-Vaidja spacetimes that reside either side of 
the brane. It was shown how the ingoing five-dimensional AdS-Vaidja spacetime was 
incomplete, and that an extension was required in order to understand some of the relevant 
cases of interest, including when the black hole mass parameter on one side of the brane 
starts from a non-zero value. The conceptual arguments presented suggested that in some 
cases graviton emission would cause the black hole mass parameter to initially decrease 
(while the brane was still inside the white hole) before increasing at later times, a fact 
that has either been ignored or has caused some confusion in the literature [41-43]. The 
case where the black hole mass was initially less than zero was discussed, and it was shown 
how here, graviton emission would always lead to the mass parameter increasing and in 
some situations becoming positive. These conceptual arguments were used to interpret the 
analytical and numerical results in the rest of the paper. 

In section 4 we solved the evolution equations in the high energy/early time and low 
energy/late time limits in order to understand the asymptotic behaviour of the solutions. 
It was shown how even though the non-Z2 symmetric parameter F decreases with time, 
the non-Z2 symmetric term in the Friedmann equation remains constant as p ^ oo, just 
like the symmetric case. Similarly, the dark radiation term C/r^, if initially large was found 
to be proportional to p. In the late time regime it was shown how both C and F rapidly 
tend to constants as the brane is cooler and much less graviton emission is occurring. 

Numerical results were then presented, showing the behaviour of each of the four terms 
in the Friedmann equation for various initial values of C and F. Graphs of C, F, 
and C~ were also given which exhibited the behaviour described in section 3.1 (and in 
figures 3 and 4), which corresponds to nonzero initial values of and C~: cases which 
have previously been ignored in the literature^. We demonstrated that in general a large 
non-Z2 symmetric term will lead to a decrease in the asymptotic ratio of dark radiation 
to normal radiation given by ej. If F is large enough then C and hence ej will become 

■^Note that in [43] they do examine one case where initially C = and F > Q, and hence C+ > and 
< 0. The brane therefore starts inside the horizon on its positive side, even though the authors have 
assumed throughout that the brane is always outside the horizon. 
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(a) Contributions to the Fricdmann 
equation where F is large enough 
to ensure that C turns negative and 
hence leads to — 0. 
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(b) Contributions to the Fricdmann 
equation where the C term turns neg- 
ative then positive. 




(c) Graph of {C/r*)/{2p){t oo) for 
a = 1.0, 0.4, 0.2, 0.02 against different 
initial values for the F term. 



(d) Graph of {C/r^)/{2p){t oo) 
against initial F term for a ~ 0.02 
and = 0.5,0.25,0.0,-0.25. Note 
the nucleosynthesis bounds. 



Figure 6: Numerical Results 
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negative which in some cases leads to i7 = at which time the universe will collapse. 
We then used the nucleosynthesis constraints on ej to restrict the possible initial values 
of the dark radiation and Y10Y1-Z2 symmetric terms, finding among other things that if the 
initial ratio > 0.1 then F must be nonzero and the brane has to be \10n-Z2 symmetric 
otherwise the nucleosynthesis bound would be violated. Interestingly, these constraints are 
all independent of M5, the five-dimensional Plank mass. 

Another interesting mechanism, besides graviton emission by the brane, that would 
affect the evolution of the non-Z2 symmetric and dark radiation terms is the possibility of 
the bulk black holes Hawking radiating. This has non-trivial effects on the evolution of C 
and F, especially when combined with graviton emissions, and will be investigated in [57]. 
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